Pseudo-proper ferroelectricity in thin films 
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We study ferroelectricity in thin films of pseudo-proper ferroelectrics such as the so-called spiral 
multiferroics. We find that this type of ferroelectricity stands better against depolarizing fields than 
conventional one. Its single-domain state can be easily preserved by metallic electrodes even in 
ultrathin films. In fact, single-domain ferroelectricty can be generated as a metastable state in the 
absence of electrodes. We also find a new regime of small thickness where unscreened films develop 
unusual multi-domain states with properties determined by non-electrostatic boundary conditions. 



Introduction.- The discovery of a new type of ferro- 
electricity caused by a cycloidal ordering of magnetic mo- 
ments [T] is one of the key events that has triggered the 
current interest in multiferroics. This type of multifer- 
roicity is realized in the rare-earth manganites i?MnO,3 
(R = Gd, Tb, Dy), and has been successfully explained 
on the basis of the inverse Dzyaloshinskii-Moriya interac- 
tion [5] . From the phenomenological point of view, these 
systems can be considered as new examples of pseudo- 
proper ferroelectrics [3] in which the electric polarization 
is (bi-)linearly coupled with the primary order parameter 
of the transition [3] . 

Much of the interest in ferroelectrics concerns their 
thin-film properties since they are at the root of a 
large number of applications (in memory devices, field- 
effect transistors, etc.). For conventional proper ferro- 
electrics this has been addressed theoretically by means 
of both first-principles calculations [5 and Landau-like 
approaches [B] . One of the main outcomes is that single- 
domain (uniform) ferroelectricity is very difficult to re- 
tain below certain sizes. Instead, there appears multi- 
domain structures which are largely the result of the de- 
polarizing fields that invebitably persist in real devices 
[7J. As regards pseudo-proper ferroelectrics, the thin film 
properties of these systems have not yet been addressed 
in the same detail. 

In this paper we show that pseudo-proper ferroelec- 
tricity proves more robust against depolarizing field ef- 
fects. We illustrate this robustness in different experi- 
mental situations. When the film is sanwiched between 
metallic short-circuited electrodes, for example, the crit- 
ical screening length necesary to keep single-domain fer- 
roelectricity is much larger than in conventional fer- 
roelectrics. In TbMn03, for example, it is expected 
~ 150A, which largely exceeds the screening length of 
many metals. The physics behind this result has noth- 
ing to do with any eventual smallness in the induced 
polarization. This can be noted in a convential ferroelec- 
tric, where no increase of the critical screening length is 
obtained by the mere reduction of its spontaneous po- 
larization. The reason is that the aforementioned lin- 
ear coupling sets a new length scale A in the problem, 



and this quantity overwhelms the typical length scale for 
the gradients of polarization in determining the above 
critical screening. On the other hand, in the absence 
of electrodes single-domain ferroelectricity can be gener- 
ated as a metastable state in striking contrast to conven- 
tional ferroelectrics. This is possible because the pseudo- 
proper mechanism does not require the vanishing of the 
polarization stiffness. Therefore, despite the depolariz- 
ing field increases this stifness, the virtual single-domain 
ferroelectric instability remains sufficiently close in tem- 
perature. In addition, the new length scale A opens a 
small-thickness regime in which non-electrostatic bound- 
ary conditions become important for multi-domain ferro- 
electricity. This is evidenced in the period of the states 
that cause the instability of the paraelectric phase in a 
film without electrodes. For film thicknesses I A the 
period decrease by reducing I as in conventional ferro- 
electrics. When I reaches A, however, it becomes compa- 
rable to the film thickness, and therefore the average out 
of the corresponding depolarizing field ceases to be effi- 
cient. A further reduction of I then produces an increase 
of such a period whose precise form depends on the non- 
electrostatic boundary conditions. This latter behavior, 
unlike in the magnetic case [5j, is rather unusual for a 
ferroelectric and, to the best of our knowledge, has been 
unnoticed so far. 

Equations of state.- The instability towards ferro- 
electricity in the pseudo-proper case can be analyzed 
from the (linearized) equations of state: 



ap - fn = -d z v, 



c x dl)r) -fP = 0, 
{e y d z y +d 2 z )V -And z P = 0. 



(la) 
(lb) 
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Here P is the electric polarization (assumed to be perpen- 
dicular to the film), V is the electrostatic potential and 
X] is the primary order parameter. The coefficient a is 
therefore the control parameter, a = a'(T — To), whereas 
A represents the bare polarization stiffness, / is the cou- 
pling constant that allows for ferroelectricity, cm and cj_ 
account for the extra stiffness of non-uniform distribu- 
tions of r), and e v is the in-plane dielectric constant. Eq. 



(lc) results from Maxwell's equations and Eqs. (la) and 
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( lb ) are derived in the appendix for the case of spiral mul- 



tiferroics [where r\ describes the transversal component of 
the magnetic cycloid: M = (0, £ cos Qy, rj sin Qy)] [9]. 

The fact that A does not vanish allows us to neglect 
derivatives of P in Eq. ( fla| ). Moreover, the essential 
physics due to depolarizing field effects can be revealed 
putting c± = 0. In this hypotetical case the order- 
parameter turns out to be insensitive to its boundary 
conditions (or, more physically, it can adjust to them 
within a zero-distance interval). This simplifies notably 
the algebra, and proves to describe correctly the behav- 
ior of the system for the case of the so-called natural 
boundary conditions d z r\ — 0. So, keeping in mind these 
restrictions, we take advantage of this model case to in- 
vestigate the point at which the paraelectric phase losses 
its stability. That is, the point at which, by decreasing a, 
there appears the first nontrivial solution of the system 
of equations ([!]) that, in addition, satisfies the remaining 
(electrostatic) boundary conditions. 

Perfect screening. - First of all, it is convenient to re- 
vise the case in which the depolarizing field is completely 
screened. In that case, the stability is lost with respect to 
the single-domain state in which both P and n are con- 
stant. This happens at the critical value of the control 
parameter a c — f 2 /A, which corresponds to the critical 
temperature T c = T + 9 where 9 = f 2 /{a'A). The in- 
verse susceptibility in the paraelectric phase then can be 



written as x e 



A 
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A- 



(T-T c )/e 



This implies a 



l l + (T-T c )/0- 

Curie- Weiss behavior in the vicinity of the phase transi- 
tion, as obtained from the first term in the expansion in 
powers of T - T c : xl 1 ~ M T ~ T c)/@- Experimentally 
this behavior is observed only in a very narrow region 
around T c (see e.g. pQ). This can be understood as a re- 
sult of a small 9 that reveals that the coupling between 
P and r] is effectively weak. In TbMnOs, for example, it 
can be estimated 9 < 1 K (T c ~ 27 K in this case) . 

Partial screening.- If the screening is not complete, 
there is a competition between the above uniform solu- 
tion of the equations of state and multi-domain structures 
that vary periodically within the film plane. Close to the 
transition point these complicated structures reduce in 
practice to a single harmonic, i.e., ~ e Thus, for 
a given k y , the functions describing the corresponding 
structure can be sought in the form n = rj^e yV cos k z z, 
P = P (r) )e tk v y cosk z z and V = V (rj )e lkyV sin k z z in 
view of the form of the equations ([I]) and the symmetry 
of the problem. Substituting in ([I]) one can see that the 
existence of such solutions implies the relation 



A + Ati- 



(2) 



between the parameters k y , k z and a (recall that cj_ = 
for the while). 

The largest values of the parameter a are obtained for 
solutions with k z <C k y . Physically this is because these 



structures produce the minimal electric field in the fer- 
roelectric, and therefore have the lowest stiffness. Con- 
sequently these structures are the most natural suspects 
of being responsible for the instability of the paraelectric 
phase as they are in conventional ferroelectrics. To move 
on we have to check whether this is actually compatible 
with the corresponding boundary conditions. 

Real electrodes.- Consider first the case of a film 
sandwiched between short-circuited electrodes. We 
model the imperfect screening in the metal with insul- 
tating dead layers of thickness d [5J and assume that 
the film thickness is I ^> d in the following. Then, to 
satisfy the electrostatic boundary conditions, the above 
solutions have to be such that 



k z l _ 
k z tan — — — £yky 



(3) 



if k y d <C 1 as we expect. The maximun a will be obtained 
for k z l <C 1, which is clearly compatible with k z <C k y 



J—h 2 

e v d h 'z- 



if d < I. We then have k% = (l + § ( Y) + • 
Substituting in ^ we can see that, below the critical 
dead-layer thickness 



3A 



Tre„ 



(4) 



where A = {Ac) 1 ' 2 /\}\, the single-domain state ap- 
pears before these multi-domain structures for the value 
a c = f 2 /(A + 4-7T j) of the control parameter. 

As we see, d c is inversely proportional to the strength 
of the linear coupling between P and r\ which, "by defi- 
nition," has to be relatively small in pseudo-proper fer- 
roelectrics. In TbMnC3, for example, the estimates 
of the spin-phonon coupling given in |10j indicate that 
A - 150A. On the other hand A, e y ~ 20 - 30 in ac- 
cordance with pp. Thus, in contrast to the conventional 
case [B], the screening of practically any electrode will 
prevent the splitting of a pseudo-proper ferroelectric into 
different domains at the transition point. 

No electrodes.- The hampering of multi-domain 
states in pseudo-proper ferroelectrics is also manifested 
for dead layers thicker than d c . Consider the extreme 
case in which there are no electrodes, i.e., there is no 
screening of the depolarizing field. In this case, the elec- 
trostatic boundary conditions are such that 



kj _ 



(5) 



for the multi-domain structures considered before. k z 
then has to be ~ tt/1 to satisfy the condition k z <C k y . 
Accordingly k y sa - ^ k j k ^ and, substituting into Eq. 
([2]), we can see that these states do not appear before 
the single-domain solution if the thickness of the film is 
smaller than A. One may then naively conclude that, be- 
low this thickness, the loss of stability takes place without 
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the formation of domains since the multi-domain states 
that normally appear in conventional ferroelectrics get 
suppressed [11] . However this is not the end of the story, 
since less conventional structures may come into play. 
We have so far ruled out structures with k z > k y be- 
cause they seem unfavorable from the point of view of 
the depolarizing field. This option, however, has to be 
reconsidered for / < A. It corresponds to solutions with 
k z l « 1, for which k y = (|^J (l + § (^) 2 + . .. J k z in 

accordance with One can see that, in fact, these so- 
lutions appear before the single-domain state for I < A, 
thus causing the instability of the paraclcctric phase. 

We note that the period of the conventional structures 
that appear for I 3> A varies ~ (A/) 1 / 2 with the film thick- 
ness (see Fig. [I]). This period becomes larger than the 
film thickness if I < A, which explains the further ten- 
dency of the system to get rid of these structures: they 
are no longer effective in reducing locally the depolar- 
izing field. Once this efficiency is lost, the factors that 
determine the subsequent behavior for I -C A are mainly 
intrinsic (i.e., the linear coupling, gradient stiffness and 
non-electrostatic bondary conditions). Thus the period 
of the structures that appear in this regime varies un- 
usually ~ A 2 jl (see Fig. [TJ. It is worth noting that 
conventional ferroelectrics only develop the former struc- 
tures whose period increases with the film thickness, as 
will be the case of pseudo-proper ferroelectrics with ex- 
tremely small A's unattainable experimentally (i.e., with 
strong couplings). 

This scenario, in which natural boundary conditions 
for r\ are implicit, holds for cj_ 7^ 0. The only changes 
are the following. The relation obtained instead of Eq. 
([2]) has two k z roots for a given k y . Accordingly the 
general solution of the linearized equations of state is 
the linear combination of the corresponding functions. 
For I > A the second k z is associated with a surface 
contribution that gives exponentially small corrections 
to the structures described before. The paraelectric in- 
stability is therefore practically insensitive to the non- 
electrostatic boundary conditions as in conventional fer- 
roelectrics. For / < A, however, the second k z becomes 
relevant to describe the z-dependence of the structure 
that appears at the transition point. But this changes 
neither the period of the structures obtained above nor 
the fact that these solutions appear before the single- 
domain state. These latter results, however, are sensitive 
to the precise form of the non-electrostatic boundary con- 
ditions. 

We also note that pseudo-proper ferroelectrics have the 
follwing remarkable property. We have seen that, with- 
out electrodes, the instability of the paraelectic phase 
implies the appearance of multi-domain ferroelectricity. 
But the single-domain state is relatively close in energy in 
spite of its depolarizing field (A + Att ~ A). Thus, by low- 
ering the temperature, it is possible to have the situation 
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FIG. 1: Period of the multi-domain structures for an un- 
screened film of a pseudo-proper ferroelectric a function of 
its thickness. A is the characteristic length associated with 
the linear coupling between the polarization and the primary 
order parameter. Inset: schematic illustration of the corre- 
sponding behavior. 



in which the free energy develops at least a local mini- 
mum about the single-domain solution. That is, a situa- 
tion in which single-domain ferroelectricity is metastable. 
This can be revealed by considering the tentative state of 
local equilibrium r/^ = — | (a — ^ , where b is the co- 



A+47T / 



efficient of the nonlinear term brf omitted so far in ( lb I , 
and computing the stiffness associated with the pertur- 
bations ?7o — > i]o + ff that could drive the system out of 
that state. It is clear that, among all possible perturba- 
tions, those associated with the multi-domain structures 
considered before are the best candidates to do this job. 
But if the control parameter is a < A + i7T (l — ^f) the re- 
sulting stiffnesses are positive. The single-domain state 
is then robust against its splitting into different domains, 
even though the corresponding depolarizing field is not 
screened (which is unthikable in a conventional ferroelec- 
tric). 

Conclusions.- We have studied the specific features 
of the instability towards ferroelectricity in thin films 
of pseudo-proper ferroelectrics. Single-domain (uniform) 
ferroelectricity can be kept in films sandwiched between 
short-circuited electrodes without requiring any excep- 
tional screening in the metal. In fact, this state can be 
generated as a metastable state even without electrodes. 
In addition, the properties of multi-domain states re- 
veal new fundamental physics related to non-electrostatic 
boundary conditions. 

We acknowledge P. Bruno, M. Civelli, E. Kats, and I. 
Paul for very fruitful discussions. 
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Appendix: Equations of state for spiral multiferroics 

In spiral multiferroics such as rare-earth manganites -RM11O3 (R = Gd, Tb, Dy), ferroelectricity is associated with 
the magnetic phase transition from a longitudinal spin-density-wave to a cycloidal distribution of magnetization. In 
the following we parametrize the magnetization as M = (0, £ cos Qy, r\ sin Qy) for the sake of simplicity. Therefore g is 
the "primary" order parameter of the transition of our interest, while the concomitant polarization P is "secondary" 
in the sense that, whithout the coupling to g, it would remain zero. The relevant coupling is due to the inhomogeneous 
magnetoelectric effect [T], which can been taken in its simplest form [2]: 

Fme = foP ■ [(M • V)M - M(V • M)], (6) 

since we are interested in distributions of polarization that vary only at relatively large scales. 

To our purposes, it suffices to consider space variations along the y and z directions and the z-component of the 
polarization only [3J. Thus, neglecting for a while the depolarizing field, the free energy of the paramagnetic phase 
can be taken as 



> =^P 2 + ^ [(dyP) 2 + (d z P) 2 ] + foP [M z (d v M y ) - M y (d y M z )} + ^M 2 + °^M 2 + *-\M 
c 



^\(d y M)\ 2 + y -\(d 2 y M)\ 2 + c ± \(d z M)\ 2 . (7) 

The appearance of the longitudinal magnetization wave M = (0, £ cos Qy, 0) can be described by assuming that 
C11 > 0. Thus, the wavevector of this modulation is Q = y/cu/g. Putting P = this structure transforms into the 
cycloid M = (0, £ cos Qy, g sin Qy) when a z + bM 2 /4 — c|/ (2g) = 2a = 0. Without clamping P this changes as follows. 
First of all, we have to consider the possibility of having M z — M z (y,z) ^ and P = P(y,z) ^ simultaneously 
Then, in accordance with free energy Q, these quantities must satisfy the constituent equations 

[A - C(d 2 y + d 2 z )] P + fo [M z (d v M y ) - M y (d y M z )] = -d z V, (8a) 
(a z + bM 2 + c l{ d 2 + 9 -d A y - 2c ± 8 2 )M z + f [2P(d y M y ) + M y {d y P)] = 0, (8b) 

where V is the electrostatic potential in the system (due to space variations of P). Putting M z — r/(y, z) sin Qy + 
8(y, z) cosQy and linearizing the above equations we get 

[A C(d 2 y + d 2 )]P- /(r? + ^{d y 9)) = -d z V, (9a) 

(a -c l] d 2 -c ± d 2 )r ] -fP = 0, (9b) 
2(a - cy^ - c ± d 2 )0 + fQ-\d v P) = 0. (9c) 

where / = /o£Q and a = [a y + bM%/4 — c 2 / (2g)}/2. Since g, 8 and P are expected to be smoother functions than 
sin Qy and cos Qy, higher harmonics have been neglected. 

We note that the function 9(y, z) describes local changes in the phase of the magnetic cycloid which, in accordance 



that ( 9c ) , are associated with the non- uniformity of the distribution of polarization. In our problem, this distribution 
is expected to vary at distances much larger than the period of the cycloid 2n/Q. Therefore, the phase of the cycloid 
remains practically unaltered and actually can be neglected for our purposes {9/g ~ k y /Q <C 1, where k y is the 



wavevector for the space variations of P). Eq. (9c I then can be omitted, and Eqs. ( |9a| ) and (9b) reduce to Eqs. (la) 
and (lb). 
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